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ABSTRACT 

A study has been made of a C-ll Instrument Flight Trainer to in- 
vestigate its capabilities as a variable stability model and simulator. 
The cockpit, with its controls and instruments, and the power supplies 
were found to be suitable for this purpose. All other components were 
not. The existing computers have been simplified to handle only lin- 
earized approximations of the aircraft stability equations of motion. 

An analog computer will have to be designed to solve the equations of 
motion in higher order differential form. In addition, systems will 
be required to provide simulated cockpit motion and visual terrain 
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CHAPTER I 



INTRODUCTION 

A variable stability model of an aircraft is a laboratory device 
which simulates the motion and response to control variation of an 
aircraft in flight. The aircraft is described by parameters such as 
mass, moments of inertia, geometry, propulsive power, etc. These in 
turn are expressed in the form of coefficients and derivatives, the 
values of which may be arbitrarily assigned and varied making the 
design of the simulated aircraft variable. When such values are mag- 
nitudes of voltage established by potentiometers, they can be intro- 
duced to an analog computer which represents, in electronic circuits, 
the Newtonian equations of motion. The analog computer solves these 
equations of motion and produces electrical signals representing the 
displacement, velocities and accelerations that the simulated air- 
craft would experience if it were in actual flight. These signals 
can be plotted on graphs, be displayed as instrument readings on a 
simulated instrument panel, provide forces on simulated cockpit 
controls, or drive a mockup model of the aircraft if it has the cor- 
responding degrees of freedom. 

If the computer inputs which represent the aerodynamic and pro- 
pulsion controls are generated by use of manual control devices in 
a cockpit mockup, a pilot can, in a sense, fly the stability model 
as if it were an actual aircraft. 

In this way an aircraft does not have to be built and flown in 
order to evaluate its flight and handling characteristics. Esti- 
mated changes in design can be made on the spot and results im- 
mediately compared. More important, the margins of stability and 
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and control can be investigated in detail without exposing personnel 
and equipment to the risk of catastrophic failure. 

An instrument flight trainer, which is normally used for the train- 
ing and maintenance of proficiency of pilots in instrument flight and 
navigation procedures, consists of a cockpit mockup with instruments 
and controls as well as navigation and communications equipment. It 
responds to the pilot's use of the controls by way of instrument read- 
ings, control forces and, in some types, motion of the mockup within 
certain limited degrees of freedom simulating actual flight. 

These are the same characteristics as those ascribed to the sta- 
bility model. The exception is that the trainer is based on the re- 
sponse of one aircraft design whose parameters and derivatives are 
fixed constants rather than any that can be variably assigned. If 
the resistances which establish these fixed parameters can be re- 
placed by potentiometers, and if the analog circuits of the trainer 
are fully compatible with the equations of motion, the trainer can 
be converted into a variable stability model and simulator. 

The subject trainer is a Link Aviation type C-ll Instrument 
Flight Trainer which has been acquired by the Department of Aero- 
nautics of the Naval Postgraduate School. The following study has 
been made to determine the suitability of this trainer for conversion 
into a variable stability model and simulator. 
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CHAPTER II 



DEVELOPMENT OF THE EQUATIONS OF MOTION 

The basic function of the flight simulator is to solve the sim- 
ulated equations of motion of the aircraft by means of analog cir- 
cuitry. This results in instrument readings indicative of position, 
attitude, and rates of change. Such results must be dimensional and 
in real time if a subject pilot is to properly respond to them. To 
provide such results, the original equations of motion must be fully 
dimensional and in real time. Since the classical analysis of air- ( 
craft stability makes use of nondimensional equations of motion in 
aerodynamic time, it is necessary to develop more suitable equations 
for a flight simulator. 

Assumptions 

The following assumptions can be made: 

a. The aircraft is assumed to be a rigid body. 

b. Earth-air coordinates are assumed to be fixed in space with 
no wind. 

c. The mass of the aircraft remains constant. 

d. The x ,z plane is a plane of symmetry. 

d. Disturbances from steady flight are small. Products and 
higher orders of velocity increments are negligible. 
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Starting with Newton's Second Law: 



Force equals the time rate of change of momentum. 

if* - 

ZFy* 



Torque equals the time rate of change of moment of momentum. 



HL 



cL hx 

IT 



= 



cLhy 

d t 




The coordinate system and sign convention adopted is shown 
in Figure 1 . 



o 




(+)X 



Figure 1 

Axes Fixed in Space and Invariant with Time 
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The components of Angular momentum are: 



din*. *[P2(2) +T>y(V)- <S}X(y)- T\jc U)]oln> 

- ? ( x 2 +Z £ ) - (5 X}fdm -ft * H 

</• hy = [<$ xfcc) - 7 yCx) + Q Z(Z) -~R y( s)J 

= -? xydm + 3 C * 2 + Z 2 ) - Ryj? 

<?t h 2 - [- P z C x ) + RxM + Ry(y) - Qz(y)~jd m 
— "? ^2 c/w — Q y 2. U m ~f 1\ (x^ ■+ '}j'~ S ) dm 

Integrating: 

hx — P Jx — ^ Jxy ~ ft J x 2. 
hy = “P Jary ■+ <fly “ft c7j2 

ha = “« J**+*I a 



15 



The sum of forces becomes: 
cLU _ 

EF Z = ») "> a z 

The sum of moments becomes: 

LL = -jf*=PI*+ Pi* - 6tJxfQj X} - RJ*,-R j*2 

Im =7F 3 =-'PJ'^-Pjzr J + QI,+'5'I 7 -K^ J - 

£N = -^jr = -P Jxi-V Jzi -QJyi- QJyi* RIz + RI Z 

Note the presence of time rate of change of moment of inertia. 

To relate the equations of motion to Eulerian axes, let axes 
X,,^ ( , and 2, be orthogonal axes fixed in the aircraft with origin 
at the center of gravity. The aircraft can then have linear and ang- 
ular velocities and accelerations, but no displacement. The veloc- 
ities and accelerations will be those which would be measured by the 
instruments in the simulated aircraft, oriented with respect to the 
three axes, while the moments and products of inertia remain un- 
changed with aircraft motion. 

Let such a set of axes have an angular velocity tO r having 
components P , 6? , and K about , y, , and 2, axes and a 
linear velocity Vy having components U , V , and W along the 
same axes. 

The angular and linear velocities are related in that Vj~ ^tXJ 5 
where f is the distance between the set of axes moving with 
the aircraft and the point, or center, fixed in space. 
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If the angular and linear velocities of the aircraft axes are 
held constant without angular and linear accelerations, the aircraft 
will still experience an acceleration. Often referred to as "cen- 
tripetal" acceleration, it is the result of the cross product U> T X 
(63 t X f )• Since V t' w t^-^, the absolute acceleration of the sim- 
ulated aircraft with constant velocity is a-t«J T XV r . 

If V T is allowed to increase in the absence of angular 

- dV r 

velocity, meaning straight flight, then CL ~'j(x • 

For the case of both linear and angular velocities with linear 



acceleration: 

where: 



(a, - “77^ +■ uy T X 



W T XV T - 



cAt 

* A 

i J k 
P 6! R 

u vvr 



= (<5V--R v)i +(RU-PW)j + (?V-QU)k 



and: 



$= ^(ulrVjtWk) = ut + y j+wk 



__ a x = U +• QW - BY 

The components of O. will be: 

a, = V - PW +■ BIT 
- \aT + p V IT 



The equations for the moment of momentum about the moving X, 3 
xf { , 2 , , axes can be expressed in the same form as for the fixed 
axes, where now the moments and products of inertia will be for the 
aircraft axes. 



where: 



h x , ='PI X - § Jx* -R Jx ? 

h y. = § ly -psjyj! 

~ -V Jic? ~Q I^ 
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Since h is a vector quantity having components h*. , , and 

, it can be treated in the same manner as the velocity vector 
, which, when rotated at an angular velocity, OJ T , gave a linear 
acceleration. However, when the moment of momentum, or angular mo- 
mentum, vector is rotated at an angular velocity the result is an 
angular acceleration. 

Thus: dB- - fiB + CD t X h 

ck t d t 

where H is the absolute moment of momentum about a point fixed in 
space, and h is the moment of momentum about the moving aircraft 
axes . 



In other words a roll must be produced by a torque; but if the 
rolling aircraft is also yawed while it rolls, additional torques 
or moments must be applied in order to maintain equilibrium. 
Expanding the cross product oo T X h : 

A £ A 

i J K 



co j X h 



P 6} R 

h *. h j, K 



- (<9h r( - R h,,) l + (Rh Z| “ PhJ j + (phy- Qh Xi ) k 

Expressing the rate of change of the absolute moment of mo- 
mentum in component form and recalling that the moments and pro- 
ducts of inertia are invariant with time: 

- RJx, ♦<$««,- 1,) -PS r‘v PRJ ^ 

- R Jyi+ PR(T X - P ! J X7 

= R Ijt P^fZj-IJ-PR Jjr P* » 5 R Jx ! 
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So far the X, , y, , Z, axis system has been arbitrary as to 
orientation in the aircraft. Now with the assumption that the 
DC, , Z, plane is a plane of symmetry, which is true for most air- 
craft, the size of the equations are greatly reduced. This is 
because there is both a positive and negative t/i of equal value 
for each X, and Z, value. 

Thus: J 'yz Jm = - 0 or Jj/z = J xy = 0 

The expanded and simplified equations of motion with respect 
to the Eulerian axes are then: 

IF* = ™[u + <S}W - rv] 

E Fy = rn [ V - PW t Kll] 

E F* = on [w + PY - $[/] 

IL = PL - R J xz + < 9 R ( I 2 - l y ) - ?®J XZ 

Em - Q Iy+ PROc-I 2 ) -R £ Jx 2 + F a Jxz 

I Ky — R 1^ — F Jxz "t F (l^ — ^ R^xz 



The preceding forces and moments are those due to gravity, 

aerodynamic and thrust forces. If equilibrium flight is defined 

as unaccelerated flight where and cO - 0 t and steady flight, 

• • . , ♦ * 

as that for which U , V , W , F , Q , n = ° , then a flight 

at a constant rate of turn is steady flight which originated 
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from equilibrium flight. 

Disturbed motion can be considered to be a change from steady 
flight, and the Eulerian axes under such conditions are considered to 
be disturbed axes. The forces can now be considered to consist of 
steady flight forces plus increments of force due to disturbances. 

The Effect of Gravity Force 

Since the force of gravity continually acts toward the center 
of the earth, its components act along the Eulerian axes which are 
not fixed in space. External moments are not affected because the 
gravity force acts through the center of gravity of the aircraft. 
However, the external forces in the three Eulerian directions are 
affected and must be determined. 

It is first necessary to find the components of the weight 
vector along the steady flight Eulerian axes, , and Z c ; 

then they are to be transferred to "disturbed" Eulerian axes which 
are rotated through Eulerian angles (p , 6 , and V' about the , 
y ,and z axes. Consider the , y 0 , and z 0 axes displaced 
from the horizontal and vertical, relative to the earth, by angles 
G 0 , (Po > and ip© . Rotation about the Z axis through t/v does 
not give a component of the weight vector about any other axis. 

From Figures 2 and 3 it can be seen that: 

OCo ~ ~ Wf 4wi Go 
Vo = W t C-O-3. Q b s&X/Y) (f 0 
Z-o = W t c^cr<L Go Cpo 
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Xc 




Figure 2 

Rotation 0 o about the y Axis 




Figure 3 

Rotation <#> about the x Axis 



These are the weight components acting along the steady or initial 
flight axes* To develop the components along the disturbed axes 
three steps are taken: First, rotation V' about the z axis, then a 
about the rotated^ axis, and then (P about the rotated x axis. 

These steps are shown in Figures 4, 5, and 6 respectively. 

x, = x co-a V + * 

y s y c^rO + X 4- 

z, = z 



Figure 4 

Rotation 'p about the Z Axis 
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X u = X, c*+9- 2,^ 9 






Z 2 “ 



2 ^ c^crl & ^ & 




■^3 - 2/* C ^ 7<> $ + z i^ vl ^ 

z 3 = z 2 c^< <P ~ 7jt *P 



Figure 6 

Rotation (p about the ^ Axis 
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Combining the three sets of equations: 

X} = X oo-j o c^o-3 ^ + y oo-3 9 -\p - H 9 

- X (oc-3 Ip l Q (p - sAasy) Ip 00*0 (p ) 

y ( tOCO 'p 0&3 (p + *4A^\ Ip xUns) 9 p) 

t Z f 00-3 9 Jlisy 7 p) 

= X ( 'V' ^A^vn & 0&<S Cp + AOcsYn p A*Uy\ cf ) 

+ y ('4^» np A04sy\ £? <O0~3 (p — 0<3 p -<£^v? (p ) 

+ H ( £> o&~3 cf) 

Substituting X* , J/ 0 , and Z* for ^ , y , and £ : 

X3 - W f 4 *' VI ^>) C -C^ & OO* Ip + ( WfCO-i #p ^0) c ^ ? "p “(w^ C-CKl Pp C^> -4^*1 6? 

2/3 “ (~ W t - 4>u>v ' &p)(oOA) lj/ -<X\ \ 9 (p — -4^Vi 3p C-4 k* cp) 

+ ( w t c-^cx? Cpo ) (&0-J 'p c-en.iL p + *dAsyi 3p (p *4 *asv\ Q ) 

■f ( Wj co-J gPp c<?</ £* p ) f c-^</ £> -4X-v> (p) 

2 ^ “ f ' " ) ( o&o ip -d^vy (9 o&~o cp ■+■ -i^V> ^ 

•+ fw^ ^ -44^ ^ *dJs*\ S cocr'J (P - o<r3 qp ^uy x (p) 

+ (w^. c#<j & 0 <-&$ q) ( C*CT3 & C*cO (f) 

The left sides of the three force equations are composed of 
aerodynamic, thrust and gravity forces. Let F* , Fy , and /> be 
the first two. The six equations of motion are now: 

LF X =■ m(U + 6}\J - KV) - ac 3 

iFy = r n ( V + RU -PW) - % 

lF a = m(w 4 PV - <$ U) “ Z 3 
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I l = pi* - k Jx, + gnarly) -?<r j X j 

£ M = d}ly +VRn*-I z )- R 2 Jxi + P*Jx 2 
I N = R I z - ? J* z + P6? (l 3 -Ir) 



• ♦ • 

The Relation Between F , £ , and R ; and , 0 , and TR 

The rotation of the Eulerian axes has been accounted for, and 
the angular velocities P , 6 } , and R have been used in the pre- 
ceding equations. These are the angular velocities about the axes 

• 4 

^3 > ^3 ♦ ar| d Z 3 respectively. The angular velocities (P , O , 
and Ip , however, are about the , and Z axes respectively 

which are not orthogonal. Therefore expressing P , $ , and R in 
terms of the latter angular velocities is complex. By resolution in 
the same manner as in the previous case of the forces, they become: 

i • 

P - Cp “ -4^1 0 

4 • 

6} = Q c-r- 4 + ny n >> cf c*ro & 

R = ^ cx-? q£/ (r — t. - • 
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The rates of change of the Eulerian angles can also be ex- 
pressed as functions of the instantaneous angular velocities by solving 
the above equations to give: 

(p — P + 6? 't-et'Ki & (p •+■ P ‘Dto- I & Cp 

6 ~ Q) C^-xj Cp - R cp 



v - «(%f) - «(£#) 



Linearization of the Equations of Motion 



The linear and angular velocities can be presented in terms of 
their initial or steady state values and disturbed increments. 
U-Uo+u F - T > o + T 

V = Vo + ^ 

W ~ Wo B = R + 'i 

Substituting these quantities into the equations of motion: 

ZF* = ™|(V + + W 0 % + <$ 0 k/- + - Vo F 0 - V/t - <?o V -'ll r) 

+• ^ ^ uw o cxkj c ^ c> &t> c^-<? 6* -*«- > i qR 

+ (f'o) -4^>i 6^ J 

Z Fl, ■= >n jjf iF + R 0 Uo + Wo/? + R 0 n- + •'*«■ - F Wo ~ W /-> - K ) 

- (ff- & 0 <Po)( c -^ a TR C4-J $ t -Wri iR £ -<2^w <^R) 

+ R @ 0 ) ( C* 6-0 "Ip 6 ? ^vi-l c^jua. ^ 

— ^ 9o <po) ( c-e-* 3 & -< 2 - <-'*•'> (f )J 

£F 2 = wjf tv- + ToVo+Vop + Je<r + T r - 6!oU v - Qo* -Kf - t«) 

+ JpAsV\ i9o) ■4WISV' \ £/ Cp ”F Cp' J 

- <L* 0 W?@ o Cpc) 'Ip -4^Vvl £ <^<2 <5^) 

- ( ^ Cpo)( c ^ <f & C^<j 
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£i- = 7 ° I* -nJxi +(($oRo -tRo? -h - i y ) 

~(Mo ■*■ ®o-p + Pof + ?%) Joct. 



Z M = % Iy + (?oRo t Ro r + Ron + v*i)( l x - 1 £ ) 

-(Ro 2 -t Z Ro~-tsi z )J xz +(? 0 l *Z?of + r*) Jjci 



L* 1 = A I g - yo J xz ■+(?,$* +#,7* -f ft? +7 ^) IjJ 

"* f^<? Ko f Ro t * fyo^X + f 0 x 2 



Disturbances from steady flight are assumed small enough that 
products and higher orders of velocity increments are negligible 
in comparison with the increments themselves. Also disturbance 
angles are sufficiently small that -<Hn is approximated by the 
angle in radians and co4. is approximated by unity. Further, prod- 
ucts of these angles may be neglected. Finally any changes in air 
density may be assumed to be zero. 
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The equations of motion can now be reduced to: 



£ Ft ~ m> 4 tfo Wo + Vc ^ + Qo ~ RoV„ - Ro s -Von) 

'*($' &(,)- c ^9 9 o (fo) Tp + f Jp C<X? <£t> <?V) 

ZFy = w [f tr -f-Ro Ue + + Ro« ~ 7^, W]> - Wo jo - Pe ur) 

- (jj* -<Wn Po) lj/ - ^ <UK> (5^ ^t^v? C-0O&O C^(p„) (fj 

Z F z = tojjw' + P®V* + 7^ it +V<> 73 - 6) 0 U 0 - $0 “■ -Hof) 

+ (fr &„)& +(jj. C*0O &o -4-W) (fo) (f ~(fr c **<9 0 c+*<f£^ 

Z l— " 'joj.yi -SI ^~(Q°Ro ■+ Pc^ +tfo/z)(l2- Iy) —(fotfc>'* tyo~P "f Rof)Jxi 
ZM-- £1^ +(FoRc +Po-p -h Pe>si)(l x -J z )-(R£-tzRo'\) J Xi 4(P 0 Z +ZPor)Jxi 
ZN - AI a ~ -fJxi -»-fPo<S?o+^or + Ro^Cly- l^)-(®oRo+ P*)Jxg 



The assumption of small disturbances not only reduces the 
equations but also linearizes them so that their solution will be 
simpler. Although this assumption seems to limit the motions to 
infinitesimal disturbances from steady flight, experience has 
shown that finite deflections can be treated in this manner with- 
out serious error. 

The equations of related angular changes have also been sim- 
plified: P =. (p - Ip & 

6? = <9 + y 

R - - <9 (p 
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and neglecting products of small disturbances: 

P = <p , & 3 R = V 

Thus the instantaneous angular velocities may be set equal to the rate 
of change of the Eulerian angles. 

The equations of motion now allow the examination of the motion of 
a body disturbed from some complicated steady flight condition such as 
combined steady state rolling, pitching, yawing motion with constant 
sideslip and forward velocity. 

In order to reduce the equations of motion still further to enable 
easier visualization, a steady flight condition is assumed with wings 
level and all components of velocity zero except U<> andW. . 

V„ To — (Qo ~ Ro = (po ~ 'Ipo — 0 

The equations of motion now become: 

IF X = m + Wo % + % i 9o + ^9 co* e 0 J 

L Fy = m[ ir + a- -V 0 - ft <4^ &o - c** (9 0 j 

L F z = rn[^ - Uo? t #,] 

HL- -j l x - A J xz 
i Iy 

Ln = A I 2 - j> Jxz 

So far the only specification as to the direction of the 
Eulerian axes has been that the y axis is a principal inertial 
axis of the aircraft. If the at axis is made to coincide with the 
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longitudinal axis of the aircraft, it is also a principal axis. This 
system is referred to as the "body" system of axes. The advantage of 
this system is that the value of Jxzis zero. However, if the sc axis 
is made to coincide with the flight direction during steady flight. 



which W 0 terms vanish. The actual choice of axes will be made later. 
Expansion of the Aerodynamic Forces and Moments 

Each of the forces x , y , and 2 and the moments L , M , and 
N can be expressed as a function of the variables upon which they 
depend by expanding the forces in a Taylor series. These series have 



where ©< , ft , and If are variables; the subscript zero indicates 
that the quantities are evaluated at the steady flight conditions; 
and all second-order terms have been neglected in accordance with 
assumption V. 

Because of the symmetry of the x , z plane, the rate of change 
of the X and 2 forces and the moment M > with respect to the dis- 
turbance velocities -p , n , and V , is equal to zero. Only deriv- 
atives with respect to U. , ur , f and their time derivatives are 
retained according to the level flight assumption. If, instead of 
this assumption, a different initial condition, such as steady side- 
slip, had been assumed, the derivatives would have to be reevaluated. 



the resulting system is called the stability or "wind axes system in 



the form: 




Of course, ) would not have been zero. 
The external forces are now expressed: 



lF x =X = 
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In accordance with the level flight assumption, the following initial 
or steady flight quantities will be eliminated when the aerodynamic 
forces and moments are substituted into the equations of motion: 
x, - mj- &o = ° 
do 

Zo + 9 0 - O 

Lq ^ o 
Mo = o 

u 0 - o 

Further, the stability axes will be used for the Eulerian axes 
system so that terms containing^ will be eliminated since there will 
be no velocity initially perpendicular to the steady flight path. 
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Substituting the aerodynamic forces and moments into the linearized 

equations of motion and dividing by the coefficients of u , v* , ^ , 

* » • 

, and /l , the equations become: 

Using the notation: 



± = Y J 

7t) dVL 



_L 21 s 

Tn d)/- ~ 




i syi _ y). 

T T 7 ~ n ,M * 

•L*) ^ U. 



J_ J 5 L_ _ 

-I* «> v 



v 



for stability derivatives. 

a = X*. m- +X{c (X i-X| ^ 4 X^ £ + Xu^ + Or 

+ Xf < 5 e ■+ JCt’ ie t X( Se 

oe oe oe 

0 + UpSl-fity ft " J ft = 3 Ar V t ^ V“ + ^ -f- ^ /i 

+ y,r* 2 ft*+ ViM + + f 3 s- £ 

0r-U o % ■+ ^9 -4*^ (9c = Z„. + Z < i % + '* r '* 

+ Z Se Se 1 - £e + Z fc ie 

- -^A = it * Li- 0 + /i +Ua + /-» 7 ° + Lj>f> 

-‘■x ' 

iL u ^ ■* L L L sL ^ + ^ ^ 

^ u. + M,x iX -t f + ^ | + + MvirtZr 

A - T + A + 7 ° + Mp f 

t M u s* * ^jL + C + ^ S«+»i‘ n U + S * 
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Here another simplifying assumption is made that the forces and 
moments acting on the aircraft are functions of the velocities and not 
the changes in the velocities relative to the air mass. This is the 
"instantaneous" view of the forces acting on an accelerating aircraft 
as if it were in steady flight with velocities at a given instant. 

This is also referred to as "quasi steady flow." 

Assuming quasi steady flow, all derivatives with respect to the 
rates of change of velocity, both linear and angular (except those 
involving Ur and ir ) and all those with respect to control deflection 
accelerations will be eliminated from the equations of motion. The 
vertical acceleration, dr , and the lateral acceleration, v , are 
retained because of the effects of downwash and si dewash, respective- 
ly, from the wing and fuselage on the tail, which may be considered 
from the quasi steady flow viewpoint. 

Consider an aircraft at two instances of time, t, , and t 2 , such 
that ^ . Assuming , the wing will experience an increase 



For small At , (//,= <*/,+ = iv; -+ dr At . 

The aircraft moves forward a distance in the time interval, 

and the tail will reach the position vacated by the wing in the time 
— - , where jl\ tail length. 



in o< at t z such that a* = 



-i 

Uo 



(U> = flight speed) 



Thus: 
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Now the change in angle of attack of the tail is proportional to 
the change in downwash which is assumed to change instantaneously in 
accordance with the quasi steady theory. 

Obviously the forces and moments on the aircraft are affected by 
the angle of attack of the tail which in turn is affected by tlr on 
the basis of the quasi steady flow method. It is for this reason that 
derivatives with respect to & are retained in the equations of 
motion. 

The effect of si dewash on the vertical tail can be treated 
similarly and therefore derivatives with respect to r are also re- 
tai ned. 

The reduced equations of motion become: 

Longitudinal : 

IA + 9o = XiA. * + t + +XV •** +Xj e Se 



- Up % + - Zm, M' + 2 ^ Sr +Z S& Se 

w + O'- + M Se S e 

Lateral : 



^ + = y v \r + +y r f+ 



• <Jx£ • 

? - -r ^ 

I X 
Jx 



= L / /+Lv-f+L n ^l+L ? f t L ^ s ~ 4 S' 



4 - jo r Vv v •+ N+. 1T+ SI + Ny, Sn 



The expansion of the determinant of the longitudinal equations 
of motion in terms of an assumed exponential root, , leads to a 
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fourth order equation of the form: 

A> V + 6 } 3 + C > 2 + D2 + £ =0 

which is known as the stability quartic. A similar quartic results 
from the expansion of the determinant of the lateral equations of 
motion. It is these equations which the stability model's analog 
system must solve. In the following section the computation methods 
used in the trainer are examined to see if they can accommodate the 
stability quartic. 

Solution by Analog Computer 

In order to simulate the motion and response of an aircraft, 
it will be necessary to employ an analog computer which will ac- 
cept the pilot's control signals as primary variable inputs and the 
stability derivatives as secondary "fixed" inputs. From the pri- 
mary and secondary inputs, simulated position, rates and accelera- 
tions of motion are computed. 

Typical schematics of the six sections of an aircraft motion 
analog computer are shown in Figures 7 through 12. 
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Figure 7 

Analog Longitudinal Motion 




Figure 8 

Analog Vertical Motion 



35 








-T 



^ r 



- <9 



-p 





Figure 10 

Analog Rolling Motion 
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CHAPTER III 



DESCRIPTION OF THE TRAINER 3 

The "Jet-Propelled Aircraft Instrument Flying Trainer," Type C- 11, 
is designed so that most of its flight, engine and communications 
characteristics are electronically computed and activated. The basic 
structure consists of a metal frame which supports the fuselage with 
cockpit mockup and the instructor's console. Components and sub- 
assemblies are mounted on sliding racks which can be pulled out from 
the framework for maintenance and adjustment. The framework itself 
is sectionalized into four main parts which are interconnected with 
cabling of sufficient length to allow the trainer to remain opera- 
tionally connected while separated for maintenance work on the in- 
teri or. 

The cockpit of the C- 11 trainer is basically a mockup of the 
F-80 jet aircraft with appropriate flight, engine, communications 
and navigation instruments, indicators and controls. The control 
stick and rudder pedals are pneumatically restricted in their move- 
ments to simulate the control loading of a jet aircraft. These 
controls are connected to pneumatic cylinders and pistons in which 
the pressure is governed by a dynamic pressure factor related to 
airspeed, air density (altitude), and Mach number. The controls 
including the throttle are also connected to potentiometers which 
provide the control commands as electrical signals. 

The instructor's console has duplicate flight and engine in- 
struments and indicator lights which monitor the function of the 
computers. There is also a ground recorder which plots the sim- 
ulated track of the aircraft over the ground calibrated to a 
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sectional aeronautical chart. Another recorder plots the glide path 
of the aircraft during simulated landing approaches. The console also 
has controls for governing the simulated communications and navigation 
environment of the trainer. 

4 

The Methods of Computation in the Trainer 

Aerodynamic and engine characteristics are computed by electronic 
relay and electro-mechanical computers. These are mounted on sliding 
racks underneath the cockpit and behind the console. 

The mathematical processes of multiplication, addition, division 
and integration are accomplished by a network of potentiometers, re- 
solvers, servos and amplifiers to produce voltages or shaft positions 
equivalent to the desired factors. Signal voltages are used to con- 
trol servo amplifiers which, in turn, control the position or speed 
of servo motors. The servo motors then control the selsyns, re- 
solvers and potentiometers used in the computation of other factors 
or directly as instrument indications. 

Longitudinal Acceleration 

Longitudinal acceleration is produced by all the aerodynamic and 
ground forces acting on the aircraft along its longitudinal axis in 
the direction of the flight path and opposite to the relative wind. 
Deceleration is produced by drag coefficients such as profile drag, 
wing flaps, dive flaps, landing gear, angle of attack and forces 
developed while in a spin. While on the ground, additional de- 
celeration of brake forces is added. These decelerative forces are 
further increased by the horizontal component of path elevation 
angle, and the factor proportional to the summation of all drag co- 
efficients is modified by dynamic pressure. The factor proportional 
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to net drag is added to the factor proportional to net thrust; and 
the result, indicating accelerative or decel erati ve force, is applied 
to the true velocity servo where its integration with respect to time 
results in true velocity. 

A voltage proportional to airborne drag is derived by the sum- 
mation of the amplitudes of voltages proportional to angle of attack, 
spin, profile drag, wing flap position, speed brakes position and 
landing gear position. These voltages, identical in phase but vary- 
ing in amplitude in proportion to their effect upon total drag, are 
added together in the longitudinal acceleration summing amplifier. 

Lateral Acceleration 

Lateral acceleration is the acceleration of the aircraft along 
the lateral axis due to all aerodynamic and ground forces except 
gravity. In the air, lateral acceleration is considered to be a 
function of dynamic pressure, engine net thrust and yaw angle. 

When the aircraft is on the ground, lateral acceleration is pro- 
portional to the difference in brake pedal deflection, rudder de- 
flection and dynamic pressure. 

In straight and level flight, the angle of yaw is zero, and 
therefore lateral acceleration is also zero. If, however, the 
aircraft is yawed, dynamic pressure and engine thrust produce 
lateral acceleration that is directly proportional to the magni- 
tude of the yaw angle. 

Computer voltage proportional to dynamic pressure, engine 
thrust, rudder deflection, differential brake deflections, and 
angle of yaw are combined. The resultant voltage, representing 
the tendency of the aircraft to skid, slip or turn is made avail- 
able to the heading, rate of turn and ball-bank computers. 
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Vertical Aerodynamic Acceleration 



Acceleration along the vertical axis with respect to the relative 
wind is produced by all the aerodynamic and ground forces in addition 
to the gravitational force. The vertical aerodynamic computer in the 
trainer calculates the lift factor of the aircraft from angle of at- 
tack, wing flap position, path elevation angle, and dynamic pressure. 
Relays and phase detectors are so arranged that the lift force com- 
ponent must equal the weight component when landing, or an enforced 
take-off will occur with resultant bumpy landing. The vertical ac- 
celeration component is routed to the pitching rate, ball-bank, 
heading, path elevation angle, angle of attack, and accelerometer 
computers and is used in their various computations. 

Normal Acceleration 

The accelerometer system uses a single pointer type indicator 
which registers instantaneous "g." A voltage from the vertical aero- 
dynamic acceleration computer, which varies in magnitude, is applied 
to the accel erometer servo amplifier which controls a servo motor. 
Coupled to the motor through a gear train is the pointer on the ac- 
celeration indicator. The magnitude of the voltage from the vertical 
acceleration computer is proportional to "g." The accelerometer 
registers between zero and plus 9.5 "g." Negative "g" is not indi- 
cated. * 

Angle of Attack 

The angle of attack is the angle between the longitudinal axis 
of the aircraft and the relative wind. Angle of attack is computed 
by adding together voltages representing vertical aerodynamic ac- 
celeration and the vertical component of bank angle, modifying this 
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sum by path elevation angle, and further modifying the whole by true 
velocity. The product of these factors is added to a voltage repre- 
senting pitch rate and is applied to the angle of attack servo ampli- 
fier, which in turn drives the angle of attack servo motor. The angle 
of attack servo functions as a velocity servo. A damping amplifier 
is included in the computer to improve the simulated longitudinal sta- 
bility. 

Angle of Bank 

Six factors are simultaneously active in computing angle of bank 
or roll. Aileron deflection and angle of yaw are the primary and most 
important factors entering into the function. Both of these are more 
effective at high air speed and are considered as varying directly in 
proportion to true velocity. The effect of unbalanced lift in a turn 
due to the difference in relative speed of the wings enters as a fac- 
tor proportional to the rate of turn. This last factor is the effect 
of a change of flight path azimuth when path elevation angle is other 
than zero. The angle of bank computer functions as a velocity servo 
system when the aircraft altitude is greater than field elevation. 

The mechanical output is used to drive resolvers in the heading, bank, 
and Zero Reader vertical pointer computers and to operate sensing 
switches in the polar cap sensor computer. 

Pitching Rate 

The pitching rate is proportional to the summation of all fac- 
tors which normally affect the pitch angle. Factors proportional 
to angle of attack, elevator deflection, wing flap deflection, speed 
brake extension, landing gear position, true velocity, and Mach 
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number are added together, and their sum is added to factors represent- 
ing net thrust and the gravitational component of vertical acceleration. 
This sum is further modified by dynamic pressure and applied to the 
pitching rate servo amplifier. The pitching rate servo functions as 
a velocity servo driving two potentiometers, one in its own computer 
system and the other which serves as a factor in computing angle of 
attack. 

Angle of Yaw 

The yaw angle is the angle between the X,Z plane of the aircraft 
and the tangent to the flight path due to its rotation about the ver- 
tical axis. Yaw is affected by rudder and/or aileron deflection. 
Voltages proportional to these factors are combined and utilized as 
the initiating voltage for the yaw servo on the heading servo chassis. 
The yaw servo drives potentiometers affecting the bank angle and lat- 
eral acceleration computers. It also provides yaw function to the 
heading computer through a mechanical differential on the heading 
servo. 

Flight Path Elevation Angle 

The flight path elevation angle is between the tangent to the 
flight path and the horizontal plane. It is computed by adding the 
vertical components of bank and pitch angles and modifying their sum 
by true velocity. This output is utilized as an initiating voltage 
for the path elevation angle servo, functioning as a velocity servo 
when the aircraft is on the ground. The servo actuates switches, 
potentiometers, and resolvers in other computers requiring the path 
elevation angle in their computations. 
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Vertical Velocity 



The vertical speed computer functions as a position servo. Its 
input is a voltage derived from the altitude computer which is propor- 
tional to the vertical component of true velocity., Its amplitude in- 
dicates the rate of vertical velocity, and its phase indicates the 
direction, either climb or dive, 

Alti tude 

The altitude system computes simulated altitude by integrating 
the vertical speed component of true velocity with respect to time, 

A servo motor, whose speed and direction of rotation are proportion- 
al to the amplitude and phase of the vertical speed voltage, gives 
a shaft position representing the vertical speed component. The 
servo motor shaft drives various potentiometers and transmitters 
which provide inputs to other systems utilizing functions of alti- 
tude in their computations. A switching system on the console panel 
is incorporated to enable the instructor to control the altitude of 
field elevation, 

Mach Number 

Mach number is the ratio of the aircraft velocity along the 
flight path to the speed of sound at a given temperature and density. 
Temperature and density are functions of altitude. The Mach number 
computer combines altitude and velocity voltage factors, making a 
comparison to a reference voltage theoretically representing the 
speed of sound, A voltage equivalent to this ratio is derived and 
applied to the pitching rate computer when the Mach number exceeds 
0 . 75 . 
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Dynamic Pressure 



The dynamic pressure is a function of air density and the square 
of the velocity. It is computed by modifying a voltage proportional 
to the ratio of air density at sea level to the air density at air- 
craft altitude by a voltage proportional to the square of the velocity. 
The dynamic pressure signal then drives a servo which positions poten- 
tiometers in other computers. The dynamic pressure servo also posi- 
tions a valve in the control loading system regulating the forces on 
the control stick and rudder pedals. 

It should be noted that air temperature or Mach number is not 
considered in this computation of dynamic pressure. This is a serious 
omission when flight velocity approaches or exceeds Mach 1. 

A sufficient number of systems have been examined for an evalua- 
tion of the trainer in Chapter IV. 
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CHAPTER IV 



EVALUATION OF THE TRAINER 

Although the cockpit is a mockup of the F- 80 aircraft, the para- 
meters and stability model in the trainer, as is, do not truly simulate 
any particular aircraft. In fact, the trainer does not really simulate 
an aircraft; it only approximates one. The primary purpose of the 
trainer is training and proficiency in instrument navigation procedures 
and not the handling of a simulated aircraft. In the interest of econ- 
omy of both cost and maintenance, the flight and engine computers are 
greatly simplified. Most functions are simply "on or off" depending 
on whether they are needed or not. Other functions are linear with 
any non-linear functions being represented by linear segments. The 
problem of solving the differential quartic stability equation, which 
is the simultaneous solution of the equations of motion, is avoided 
by this simplified "curve following" technique in the trainer. 

Most of the hardware in the trainer consists of communications 
and navigation equipment which will not contribute to the variable 
stability model and can be discarded. 

The cockpit mockup with its controls and instruments can be 
the basis for the simulator portion of the variable stability model. 
Only versatility modifications are necessary to make it easier for 
different types of pilot controls and indicator presentations to be 
evaluated in future projects. 

It will be necessary to design and construct an analog computer 
which will be able to solve the equations of motion that have been 
developed here. Although none of the existing computers can be 
utilized, they will furnish a supply of component parts and sub- 
assemblies which might be used in constructing the analog circuits. 
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The power supplies presently in the trainer appear to be adequate 
to support any new equipment which will have to be constructed. They 
accept either single phase, 230 volt, 60 hertz or three phase, 208 
volt, 60 hertz and produce a comprehensive assortment of AC and DC 
voltages. 

The parameters which are fixed in the existing computer systems 
are buried in the circuits and often in several different places due 
to the piecemeal methods of computation. This precludes bringing 
them out to a panel and replacing them with only one variable poten- 
tiometer each. This means that a potentiometer panel to handle the 
variable parameter inputs can only be constructed in addition to the 
new analog system. 

The pneumatic system which provides control force reactions to 
the pilot is controlled by a factor which is directly proportional 
to the computed dynamic pressure. In order to provide more accurate 
control force simulation, this factor should be related to the re- 
spective hinge moments that are developed and thereby subject to any 
parameter variation. This is not considered to be too important 
since it is anticipated that most of the future projects using the 
stability model as a simulator will be concerned with full power 
control systems. 

The trainer also lacks two more systems which are essential in 

an effective simulator; a simulated visual terrain reference out- 

5 

side the cockpit and a cockpit motion system. 

The visual reference is absent because the trainer is concerned 
ed. strictly with instrument flight. Normally, however, a pilot re- 
lies to a great extent on the visual cues he receives from the 
geography over which he is flying as to attitude and direction. 
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This can be simulated by a display on a projection screen which can be 
a television monitor. The p^c+u^e displayed can be a signal from a 
television camera mounted with six degrees of freedom over a scale ter- 
rain model while the relative motion between the two is governed by the 
analog solution of the aircraft equations of motion. By another method, 
the picture displayed could be generated by a graphical spatial analy- 
sis computer which can provide the illusion of three dimensional motion 
over a grid surface. What this latter method lacks in realism it makes 
up for by having infinite geographical limits. 

To be accurate, the cockpit of a simulator must also be provided 
with the same six degrees of freedom as the actual aircraft regarding 
accelerations; three angular and three translational. However, this 
is precluded by space and cost requi rements short of the real aircraft. 
The economy of the trainer is such that the^e is no cockpit motion at 
all. This omission leaves dependence on visual cues only and removes 
the pilot's tactile cues of acceleration which lead the visual cues 
by a phase difference factor. The visual cues lag because they de- 
pend on displacement, the second integral of acceleration. The tac- 
tile cue is the actual acceleration. 

There remains the spatial limitations imposed upon the excur- 
sions of the simulator cockpit. A dependence on visual cues only 
requires considerable displacements even when a scale mode 1 is used. 

Because the tactile cues lead, the pilot uses them as a sort 
of early warning for the impending visual changes. Continued motion 
is not necessary because after the initial tactile cue the visual 
cues take over. Such as been established in tests of various sim- 
ulation methods at the Ames Research Center, Moffett Field, Calif- 
ornia. It was further found that instead of a true acceleration 
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cue an<^ exponential washout or decay of such a cue will provide a suf- 
ficient tactile cue* The actual duration of the acceleration can be 
reduced. The "lurch" felt by the pilot becomes a "nudge" as the ex- 
ponential time constant is reduced, and the actual excursions required 
by the simulator cockpit can be kept within reasonable proportions. 

When necessary, certain prolonged accelerations, such as the 
take-off roll, can be simulated by tilting the simulator cockpit 
slightly while keeping the relative visual attitude cues as if still 
level . 

The trainer can be modified to incorporate a tactile cue system 
by mounting the simulator cockpit frame on a system of hydraulic jacks 
whose valves are actuated by signals from the analog motion computer. 

A general schematic of the final form of the variable stability 
model and simulator is shown in Figure 13. 
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Figure 13 



Aircraft Variable Stability Model and Simulator 
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CHAPTER V 



CONCLUSIONS 

The C-ll Instrument Flight Trainer cannot easily be converted into 
a variable aircraft stability model. However, the Trainer does have a 
suitable cockpit mockup with controls and indicators which can serve as 
the simulator portion of the stability model. The existing computers 
in the trainer will have to be replaced by an analog system which is 
capable of solving the stability quartic of the aircraft dynamic equa- 
tions of motion. The detailed design of such an analog is recommended 
as a subject for another thesis as a sequel to this study. 

The removal of the navigation, communications and computation 
equipment which is not needed for the stability model will make room 
for some of the new systems. The additional external features which 
will have to be constructed are: a potentiometer panel, to provide 
the inputs of the aircraft parameters which make up the stability 
derivatives necessary to the equations of motion; recording equipment, 
to plot the time histories of the resulting factors of motion; a cock- 
pit motion system, to provide tactile cues to the pilot; and an out- 
side cockpit display, to provide the pilot with a visual geographical 
reference. 

The existing power supplies are vacuum tube types which could 
possibly be utilized by the new equipment unless solid state circuits' 
are preferred. Solid state systems would be more reliable and re- 
quire less cooling than vacuum tube systems. 

The variable stability model and simulator will be a very useful 
tool for the investigation of stability and control problems. Transfer 
functions for human pilots as well as automatic controls can be evaluated. 
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APPENDIX I 



SYMBOLS USED IN C-ll INSTRUMENT FLIGHT TRAINER SCHEMATICS 

B Dynamic pressure substitute. 

M Mach number. 

V p Airplane velocity along the flight path. 

Angle of attack. 

Fraction of wing flap deflection. 

Si.fr Fraction of full wheels down position. 

£^1 Fraction of dive flap deflection. 

^e/, Fraction of full right toe brake. 

Sbji Fraction of full left toe brake. 

Tw Net thrust. 

^ Total longitudinal aerodynamic and ground force. 

dz. Total longitudinal acceleration. 

Y Angle of yaw. 

Rudder deflection. 

fy Total lateral aerodynamic and ground force. 

Total lateral acceleration. 

Fl<r Vertical ground force on aircraft. 

F z Total vertical aerodynamic and ground force. 

Total vertical acceleration. 

S e Elevator deflection. 

a. 2er Acceleration due to ground forces. 

Pitching rate. 

Aerodynamic pitching moment. 

Iy Aircraft moment of inertia about the lateral axix. 
ft Acceleration due to gravity. 
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Flight path elevation angle. 
<P Angle of bank, 
h Altitude. 



Constants used in C-ll Instrument Flight Trainer: 



K* 



Ky 



k. 

k 3 

Kv 

K f 

Kc, 

Ktc 



-7.8 for between -10 and 13.8 degrees. 
-26 for between 13.8 and 20 degrees. 
-8.8 for between -10 and 13.8 degrees. 
242. for between 13.8 and 20 degrees 
18.2 
45. 



2.1 

32. 

750. 

Adjustable in the range 1.0 to 0.25 
718. 

958. 

323 



0.4 



K|, = 57.3 
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APPENDIX II 



C-ll INSTRUMENT FLIGHT TRAINER COMPUTATION SCHEMATICS 



ANGLE OF ATTACK COMPUTER 
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LONGITUDINAL ACCELERATION 
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